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Abstract 



Junctions on cosmic string loops give rise to the proliferation of sharp kinks. We study 
the effect of this proliferation on the gravitational wave (GW) signals emitted from string 
networks with junctions, assuming a scaling solution. We calculate the rate of occurrence 
and the distribution in amplitude of the GW bursts emitted at cusps and kinks in the 
frequency bands of LIGO and LISA as a function of the string tension, the number of 
sharp kinks on loops with junctions and the fraction of loops in the cosmological network 
which have junctions. Combining our results with current observational constraints, we 
r~| \ find that pulsar data rule out a significant number of kinks on loops for strings with 

tensions Gfi > 10~ 12 . By contrast, for smaller tensions current observations allow for a 
large number of kinks on loops. If this is the case, the incoherent superposition of small 
bursts emitted at kink-kink encounters leads to an enhanced GW background that hides 
■ the strong individual bursts from kinks and cusps. 

o\ '■ 1 Introduction 
O 

The gravitational waves (GW) emitted from cosmic string and superstring networks may enable 
the detection of strings with very low tensions using the LISA space interferometer [Tj. In this 
paper, following up on previous work [21 E], we study the effect of the presence of 3-way junctions 
(or Y-junctions) on the GW burst signal emitted by a cosmological network of cosmic (super-) 
strings. Such junctions are thought to be a generic feature of cosmic superstring networks, 
which contain F- and D- strings as well as bound states made up of both types. Junctions 
can also exist, in cosmic string networks formed in symmetry breaking phase transitions - the 
simplest examples being Abelian-Higgs strings deep in the type I regime or Z3 strings. 

In [3] we have recently studied the dynamics of kinks interacting with Y-junctions, and 
found that the presence of junctions on a closed loop dramatically increases the number of large 
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amplitudqj kinks on the loop. That is, kinks proliferate. Indeed, when a kink reaches a junction 
on one of the strings it produces three daughter kinks — a reflected kink and two transmitted 
kinks. The daughter kinks often have an amplitude (or sharpness) that is comparable to that 
of the incoming kink, leading to an exponentially increasing number of large amplitude kinks 
on the loop. As we argued in [3], by the time proliferation ends the number of large amplitude 
kinks on a loop is generally several orders of magnitude larger than the number of sharp kinks 
expected on loops without junctions. 

Here we study the effect of this proliferation on the GW burst signal emitted by string 
networks with junctions, focusing on models in which the different strings in the network all 
have similar tensions fi. Bursts from cusps and kinks on strings give rise to high frequency GW 
signals that are superimposed on the low frequency spectrum of the string, and the waveforms 
of GW bursts on loops with junctions have been calculated in [21 SJ - In this paper, we calculate 
the rate of occurrence and the distribution in amplitude of the GW bursts emitted at cusps and 
kinks in the frequency bands of LIGO and LISA, as a function of the string tension Gfi, the 
number k' of sharp kinks on loops with junctions, the fraction q of loops that have junctions and 
the intercommutation probability p. To do so we will use a simple 'scaling' model for the cosmic 
string network. We also calculate the stochastic background generated by the superposition of 
overlapping bursts. 

On a given loop, an individual burst emitted by a cusp is stronger than a kink burst. 
However, to determine the contributions of cusps and kinks to the high frequency GW burst 
signal in a cosmological network of cosmic strings, a more detailed analysis is required. This 
was done by Damour and Vilenkin in [I] for 'standard' loops without junctions, where it was 
shown that bursts from cusps provide the dominant contribution to the GW signal. Here we 
generalize their analysis to take in account the effect of junctions. 

The paper is organized as follows. In Section [2J we describe our scaling model for the network 
evolution of loops containing junctions. We compute the amplitude of bursts emitted by such 
a network in Section [3J and compare this with the sensitivity of LIGO and LISA. We then 
turn to the calculation of the stochastic background created by the superposition of bursts in 
Section HI which is also affected by the kink proliferation. Finally, in Section we combine our 
results from Sections [3J and H] with current observational constraints and discuss the prospects 
for observation in the different regimes of parameter space. 



2 Network evolution 

The GW signatures from cosmic string networks depend sensitively on the details of the network 
evolution, and particularly on the loop distribution. Predictability is therefore still limited since, 
even in the case of standard cosmic strings with no junctions there is still significant uncertainty 
and disagreement on, for example, the loop distribution on the smallest scales (see, however, 
[S]). However, for such networks there is broad agreement that on large scales the network 
reaches a scaling solution, though its dependence on the intercommutation probability p of the 
strings is still not fully understood. 

In this paper we will follow closely the approach developed by Damour and Vilenkin in pi] for 
networks of loops without junctions. Hence we begin this section by recalling the assumptions 
made there before developing a generalization applicable to networks with junctions. 

lr rhe amplitude is defined as follows Working in Minkowski space and in the standard conformal temporal 
gauge, the position of a string is given by x(cr, t) = | (a(<r + t) + h{a — t)) where a' 2 = b' 2 = 1. The amplitude 
of a discontinuity (or kink) in, say, b' is defined by A[b'] = i 1 1 1> /+ — b'~||. 
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2.1 Scaling solution 

Central to the GW calculations for a network of cosmic strings without junctions (e.g. [UE1E1IE]), 
is the fact that the network reaches a scaling solution in which energy is lost from the network 
by loop formation. At time t the loops are all taken [I] to be formed with size 

L = at, (1) 

that is the loop production function is taken to be a ^-function (see, however, P, E3 EI] ) • There 
is no agreement on the value of a, but we shall follow [I] and adopt the simplest scenario in 
which it is set by a crude estimate of the efficiency of gravitational backreaction for the damping 
of the small scale wiggles; 

a = TGfj. , where T = 50 . (2) 

In what follows, we keep both parameters a and Gfi explicit in our formulae so that their 
origins can be traced more easily. However, when plotting our results, we will use ([2]) to express 
everything in terms of Gfi. 

Once formed, the loops decay by gravitational radiation meaning that at any given time t 
there is a distribution of loop lengths L, characterised by a number density n(L,t). Following 
[I] we also approximate the loop length distribution by a 5-function at time t (see, however, [7] 
that goes beyond this approximation and shows it does not affect order of magnitude estimates) 
so that the number density of loops of size L = at at time t is given by 

n(t) =p~ 1 a~ 1 r 3 , (3) 

where p < 1 is the reconnection probability. Notice that the typical distance between loops is 
approximately n _1//3 = p 1 l' i a~ 2 ^L ^> L. 

We now extend this discussion to networks with junctions. The evolution of such networks 
has been studied in the recent literature using both numerical and analytical approaches P21 [THJ, 
EEU E3, E3 EE] • In [T2] (see also [To])) the 1-scale model of cosmic string evolution developed in 
[T8 | [Till |2"U] has been extended to string networks with junctions, and the resulting 'modified' 1- 
scale model equations generically have scaling solutions. One therefore expects the network with 
junctions to scale, and this has indeed been observed in numerical simulations [US, El EE fTTj . 

In this paper we focus on networks of strings with junctions in which the different strings 
in the network have similar tensions /i. In that case, in the scaling solution the number density 
of loops without junctions is as given in Q, 

n{t) = p^a'H- 3 . (4) 

Here p is again the intercommutation probability which we assume is the same for all the strings. 
On top of these loops, however, there will be a set of loops with junctions (not considered in 
[T2l [To] but seen in simulations [21] )• Consistent with the scaling solution and as we argue 
below, we will take the number density of such loops n'(t) to be a small fraction, q < 1, of the 
number density of loops without junctions; 

n'{t) = q ■ n{t) = q (p" 1 ^"^" 3 ) . (5) 

2.2 Formation of loops with junctions 

To understand the formation of loops with junctions, and in particular the typical size of such 
loops, consider a network with three different strings which we differentiate by giving them a 
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colour; say red, black and blue. A junction is a point at which these three different strings 
meet. Suppose that a string of one type (blue on Figure d]) collides with a string of another 
type (black). Depending on the underlying theory and on the conditions under which this 
happens (velocities, angle,...) [22], these two strings may become linked by a third string (red), 
thus generating two junctions. If the size of the red string remains much smaller than at, the 
characteristic size of the smallest wiggles, it is expected to shrink again to zero. By contrast, if 
the length of the red string becomes comparable to at, it becomes possible that both the blue 
and the black pairs of segments intersect again and form by exchange of partners a loop with 
junctions at the scale at. □ Therefore, we assume in the following that all loops with junctions 
are produced at the same size L = at as loops without junctions. 




Figure 1: Formation of loops with junctions 

Since such loops can only form around junctions whereas standard loops can form anywhere 
along a string, the number of loops with junctions must be smaller than the number of those 
without junctions. In other words, we expect q defined in Eq. <^ to be smaller than 1. Fur- 
thermore, in the scaling regime where t sets the only scale in the problem, we expect q to be 
time-independent. 

3 Gravitational wave bursts 

In the previous section we have described the formation and evolution of a population of loops 
with and without junctions in a network of cosmic strings. In this section we determine the 
GW burst signal emitted from such a network and study its detectability in the LIGO and 
LISA frequency bands. 

This amounts to a generalization of the analysis of [I] (to which the reader is referred for 
more details) that takes in account the fact that some of loops in the network contain a very 

2 if the dynamics allows for many such red links to grow much beyond the typical size at, one expects that 
this will lead to the formation of a highly interconnected net, a possibility that we have discarded. 
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large number of kinks. The different steps in this analysis are the following. One first specifies 
the number of kinks and cusps on an individual loop of size L, both with and without junctions. 
Then one determines the rate of GW bursts dN(z) generated by cusps and kinks from loops 
between redshifts z and z + dz. The integral of this over redshift up to z m gives the total rate 
of burst events N(z m ) of different kinds produced by loops as a function of z m . By inverting 
this relation and substituting this into (El) one obtains the amplitude of bursts as a function of 
the rate. Finally we compare this with the sensitivity of the different G W interferometers and 
discuss to what extent these experiments will enable one to place constraints on the network 
parameters for the specific network model adopted here. The stochastic background generated 
by overlapping bursts will be discussed in section HI 

3.1 Waveform 

For a burst emitted from a loop of typical size L at redshift z, the logarithmic Fourier transform 
of the observed amplitude of individual GW bursts at frequency / is given by the following 
order of magnitude estimate, 

h ^~wrh?W) ^-m/.m), to 

where to is the age of the universe and the angle 6 m is given by 

9 m (f,L,z) = (f(l + z)L)- 1 ' 3 . (7) 

The Heaviside function enforces that the burst waveforms are only valid in the high frequency 
domain 9 m < 1, and finally the value of is given by 

• = 1/3 for a burst emitted by a cusp. 

• = 2/3 for a burst emitted by a kink propagating on one of the strings. 

• = 1 for a burst produced at a kink-kink encounter or when a kink crosses a junction. 

One should note that in all these three cases, numerical factors of order 1 arising in the saddle 
point calculation of the waveform (see e.g. [I]) have been omitted in For kinks, there is 
an additional factor equal to the amplitude of the kink itself, defined above. In Appendix lAl 
however, we show that only kinks of amplitude 0(1) contribute to the GW signal and hence we 
use (jSD in the remainder of the paper. One should also note that the saddle point calculation 
of the burst waveform yields expressions containing the second derivative of the left- and right- 
moving waves along the string, a" and b", which in Eqs. and (j7|) have been approximated 
by 

a", b" = C(l/L) (8) 

for a loop of length L, as in [H [?l E]- In this paper we consider loops with a large number k' 
of kinks and hence one can ask if this approximation is still valid since there is now a second 
length-scale L/k' in the problem. Despite that, we leave Eq. (jHJ) unchanged (see, however, [U]) 
since we see no particular reason why the curvature of the wiggles between two consecutive kinks 
should differ from the usually assumed 1/L curvature of kinkless loops. Indeed it is straight- 
forward to construct loops of length L with numerous kinks but arbitrary curvature in between 
the kinks, showing that the number of kinks and the curvature scale are not in principle directly 
related. We therefore use Eq. (jHJ) and the corresponding expressions Eq. (jHJ) and (JZJ). % 



3 Ultimately, however, numerical simulations of loops with many kinks must be carried out in order to 
determine the typical order of magnitude of a", b". 
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3.2 Number of cusps and kinks per loop 

Consider first the population of "standard" periodic loops with no junctions. Let c (resp. k) 
denote the average number of cusps (resp. kinks) per oscillation period T L = L/2 of a loop of 
length L. Then the average number of cusps or kinks per unit spacetime volume on such a loop 
is given by 

c 

^cusps - T n ( t ) > ( 9 ) 

k 

V kinks ?t -n(t) (10) 

where, as in the following, we drop factors of 2 since we are interested in order of magnitude 
estimates. On loops without junctions, one typically expects c ~ 1 and k ~ 1. 

The situation is more complicated for the remaining loops with junctions since these do not 
evolve periodically and, as we have shown in [3J, the number of kinks on loops with junctions 
rapidly proliferates. This proliferation is limited, either by junction collisions, gravitational 
back reaction or because of the inner structure of the string. 

In this paper, we assume that the typical timescale T pro uf of kink proliferation is much 
shorter than the typical lifetime of the loop Ti oop . Moreover, we assume that the latter is not 
affected by the presence of many kinks, that is to say that most of the power radiated by a 
loop comes from the low frequency modes so that the total power is essentially independent of 
the number of kinkfl Finally, we also assume that the typical timescale T roun d ff over which 
gravitational backreaction rounds off kinks is larger than Ti oop so that a kink propagates and 
radiates during the whole life of the loop. All these assumptions can be summarized as follows 

Tprolif 7~loop Tround off ( 1 1 ) 

The first inequality of (fTTj) is for instance realized if junction collisions end the proliferation 
process. In this case one of the strings in the loop unzips when the junctions collide, and we 
assume that as a result of this the loop splits into two separate loops [23] • These loops have no 
junctions and hence evolve periodically. However, each of these contains a large (and constant) 
number of kinks k' ^> 1 propagating on the loop. The uncertainties on the exact value of k! 
make it natural to consider it as a free parameter in our analysis. By analogy with the case of 
standard loops, we define d as the number of cusps per interval of time L: 

^-n(t) (12) 

^n(t). (13) 

Since the presence of numerous kinks is expected to inhibit the presence of cusps [21] , we expect 

d < c (~ 1) , (14) 

though our results below hold provided c'<c. 

We have estimated the order of magnitude of k! in the scenario in which kink proliferation is 
limited by junction collision by numerically evolving loops with three strings and two junctions 
[3J. These simulations show that, after a typical time T pro uf ~ 10L, junctions collide, and at 
this time the total number of kinks on the loop is ~ 2> Tvrolif ' L . However, not all kinks contribute 

4 The validity of this assumption when the number of kink is large will be investigated in future work. 
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to the signal given in Eq. as we show in Appendix A, the signal is dominated by large 
amplitude or 'sharp' kinks. Taking k! to be the number of sharp kinks, defined as kinks with 
amplitude A>0.25, it was found [3] that each string that makes up the loop contains of order 
10 3 kinks moving in one direction (left or right). This number must be multiplied by 3 - the 
number of strings in the loop - and by 2, to account for the left moving and the right moving 
kinks, so that kl ~ 10 4 . In cases where junctions collide later, say when T pro uf ~ 15L, one 
easily reaches k' ~ 10 6 . In fact, for more complex loops with junctions, the amplification of 
the number of kinks might continue long after the first junction collision e.g. because of the 
formation of new junctions [221, so that values of k! ^> 10 6 do not seem unreasonable. 
In the following we concentrate on models for which the combination qk' satisfies 

qk' > 1 . (15) 

We believe this is a realistic class of models, since q < 1 whereas one always expects k' ^> 1. As 
can be seen from Eq. (1131) . from the point of view of the rate of kink events, the subnetwork of 
loops with junctions behaves as a network of standard loops with an effective number of kinks 
qk'. Hence in the class of models we consider the average number of cusp events is dominated 
by loops without junctions, since q < 1 and c 3> c', whereas the dominant contribution to kink 
events comes from loops with junctions. 

3.3 Amplitude of cusp and kink bursts 

From now on all quantities will be expressed in terms of redshift z. For simplicity we also use 
the same interpolating functional between t and z as in [I] . Thus 

t~t <fi(z) where tp^z) = (1 + z)~ 3/2 (l + z/z^ 1 ' 2 (16) 

where z eq = 10 3 ' 9 is the redshift at radiation/matter equality, and to = 10 17,5 s is the age of the 
universe. Below we also need dV(z), the proper spatial volume^ between redshifts z and z + dz; 

dV{z) ~ 10 2 t 3 2 2 (l + 2)" 13/2 (1 + z/z eq )~ 1/2 dz . (17) 

Let rfA r cusp /kink(2 ; ) be the rate of bursts generated by cusps or kinks that are emitted between 
redshifts z and z + dz and reach us today. Then 

dN cusp (f,z) = ^(l + z y l u cusps (z)dV(z) (18) 
dNuAf,z) = e m (l + z)- l u'^(z)dV(z) (19) 

where 9 m {f,z) is defined in Eq. ([7]). The first factor in each equation is the beaming fraction 
of the cusp/kink, namely the probability that we are inside the set of directions of emission of 
a burst produced at redshift z. The second factor is the usual time dilatation factor. The third 
factor, the number of cusps/kinks per unit spacetime volume, has been estimated in Section 
13.21 and is given by 

^cusps(^) ^ c (p^crto V/(^)~ 4 ) , ^kinks(^) - Qk' (p^a-X Vi(^)" 4 ) • (20) 

5 In [7], it was shown that using these rather than more accurate interpolating functions (taking into account 
the recent change from matter to vacuum energy domination) does not affect order of magnitude estimates. 
6 The numerical factor of 10 2 in (|17|) approximates an exact numerical factor that can be found in [J. 
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For each cusp or kink signal one can solve for the smallest redshift (2™ p or z^ nk ) needed to 
ensure a certain detection rate (N cusp or iVkink)- One should bear in mind that events occuring 
at a rate smaller than 1 yr" 1 are not likely to be observed by experiments. Dropping the 
kink/cusp label, we therefore need to solve 



N 



dN(f,z) 



(21) 



This integral is dominated by the largest redshift so it follows directly from eqs. (I18p - (l20p that 



TV ~ 

1 v cusp — 



io 2 a-^iftoy^itopy'c (c sp ) 3 (i + c sp r 7/6 (i + z c r/ze q ) 



vll/6 



10 2 a-Wtftor^M-W (^ nk ) 3 (l + e k )" 4/3 (l + 



.5/3 



(22) 
(23) 



'kink — \j »vj i-" \~ m ) \^ ' ~m / r 1 ~m 

Note that in each case N increases monotonically with z m . Inversion of the above gives 
z m = Zm(f, N) so that, using eq. (ED, the observed signal is 



h(f,N) = h(f,z m {f,N) 



(24) 



The amplitude h decreases with z m and thus with N. Its explicit form reads for cusps and 
kinks respectively, 



h{f,N, 



cusp. 



GfM {10' 2 N cusp t 
Gn (lO- 2 iV cusp t 
Gfi (lO~ 2 iV cusp t 



-1/3 



a 



-2/9 



(fk 



-5/9 1 



-8/11 
-5/11 



a 



-14/11 



(fk 



c/p) 

-9/11, 



1/3 



X 



cusp 
m 



z/pfl 



11 



, Zeq » ^ > i; 



a- 6/11 (/*o)- 7/U (c/p) 5/11 ^ 1/a (« 



cusp 
P > 2 



M/^kink) ~ < 



-1/3 
-7/10 
-2/5 



(25) 



a- 4/9 (/to)- 7/9 (#'/p) 



G/i (lO~ 2 iV ki nkto 
G/i (lO" 2 iV kink t 



7 / 9 (qk'/ P y/ 3 
13/10 (fk)- 9 / 10 (qk'/py/ 10 



v kink 



«1) 



(z eg » z^ k » 1) 



(:'«/ ^ 



(26) 

where we have dropped the Heaviside function in eq. since it will play no role in the frequency 
bands and string tensions of interest. 

We now argue that this by now standard procedure should be modified by the introduction 
of a lower redshift cutoff z c . This arises from the observation that the above analysis relies on 
a continuous description of the population of loops in the network (through the density n(z)) 
which, in the context of the model described in Section El is only valid above a certain scale 
(roughly n(z)~ 1/s ). Therefore, this approach does not apply to redshifts so small that there are 
not enough loops in a sphere of the corresponding radius to allow for a statistical approach. 
Hence ;z£ usp and z kmk can be estimated by imposing 



N(z c c 



cusp \ 



n(z)dV(z) 



N(z, 



kink\ 



Q 



n(z)dV(z) 



1 . 



(27) 



The value of z c only depends on a, p and q, but not on the number of cusps or kinks per 
loop: indeed, from eqs. (I27p . and since z c <C 1, it is straightforward to show that 



cusp 



/ pa\ 

uov 



1/3 



kink ^ 
— 



[lQiq, 



1/3 



(28) 
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Note that since p < 1 and a < 10~ 5 , one indeed has z^. nsp <C 1. From z kmk we must impose 
pa <C q for consistency. We also emphasize that z c intrinsically depends on the density of loops 
in the network and not on the type of burst. The difference between cusps and kinks in ( 1251) 
is a consequence of the fact that one expects to observe cusps from the standard subnetwork 
(with q — 1) and kinks from the network of loops with junctions (q < 1). 

Hence if for one or both types of burst event z m <C z c , then the probability of there being 
one cosmic string loop inside the spheral is <C 1. In this case, no matter how large the number 
of cusp/kinks per loop, the probability of observing anything at all is small. For most values of 
the parameters and of the rates considered in [I] , the corresponding redshifts z m are above the 
cutoff value z c so their results are unaffected by this discussion. However, when the number 
of cusps or kinks per loop is very large, then in part of the range of rates of experimental 
interest one finds z m < z c (because ZfYi IS cl decreasing function of k' as can be seen from (I23| ). 
Therefore, instead of ([Ml), one should write 

h(f, N) = h (/, z m (f, AO) 6 (z m {f, N) - z c ) (29) 

when comparing the predictions of our models with observation. 

For fixed parameters (a,p,q,c,k f ), the lower redshift cutoff z c turns into a lower cutoff on 
the rate of events and an upper cutoff on the amplitude. Indeed, on decreasing the size of the 
sphere around oneself, the rate of events decreases and the amplitude of the signal increases. 
However, if at some point the sphere contains no loops then reducing its size further obviously 
has no effect - the GW signal vanishes. Specifically, from eqs. (j22|) . (123"1) together with 02 8|) the 
lower cutoff on the rates is given by 

« - k ' G*0 ((dps) ■ < 31 > 

Note that these are independent of q and p. However, one sees that kink proliferation increases 
the minimal rate. The corresponding (k'~ independent) cutoff on the amplitude h(f, z c ) can be 
obtained from (JBJ) and is given by 

*S5 = 102/3 ($) ((i^vs) « ,/3 • < 32 > 

*SS = 10 2/3 (||) ((i^Ts)*" 3 - (33) 
In the parameter range that we will consider, /i™sp ^ ^Snk- 

3.4 Detectability of Bursts 
3.4.1 Sensitivity of the instruments 

Whether or not a burst is detectable depends on the sensitivity of the experiment. We now 
estimate the optimal signal to noise ratio (SNR) of the bursts following [1]. The true SNR p is 
defined as 

2 _ f +co df houtifthl^if) 



, = k _, ^ Jr - (34) 



7 Note that the loops that we consider are basically pointlike on cosmological scales. 
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where h temp (f) is the (logarithmic) Fourier transform (FT) of the best template, h out (f) the 
(logarithmic) FT response of the detector and S n (f) is the (one-sided) spectral noise density. 
Under ideal circumstances — that is, when the detector output equals the actual signal arriving 
at the detector and when the template matches exactly this signal — eq. (j34p reduces to 

" = k p^i) ' (35) 

Going beyond this idealized situation requires a much more detailed analysis such as those 
performed in , but we limit ourselves to an estimate based on f[3"5l) which suffices to give 
an idea of the detectability of the effect of kink proliferation on the GW burst signal, at least 
when p ^> 1 and p <C 1. One should bear in mind however that values p ~ 0(1) in our analysis 
do not allow one to draw conclusions about the detectability of the signal. 

Using a signal waveform h(f) = B\f\~^e 2m ^, eq. 0351) can be rewritten as p 2 = \B\ 2 J +o ° ^s(f) 
where s(f) depends on the type of burst (the value of /3) and on the instrument. Physically, the 
value f c at which s(f) reaches its maximum is the frequency around which the instrument is 
most sensitive to that type of burst. We can use this frequency to replace \B\ in the expression 
of p 2 by h(f c )f~P to find an expression for the signal to noise ratio p = h{f^)f~^{jQ°° / s (/)) 1/ ' 2 - 
This can be written as 

p = ^t- ( 36 ) 

Note that even though the SNR is conveniently expressed in terms of the signal at a single 
frequency f c , its computation is really a noise weighted integral over the entire frequency band 
of the instrument. 

Since f c and do not depend much on the kind of burst, using one set of values per 
instrument only introduces an error of a few percent on the SNR that will not affect our or- 
der of magnitude estimates. For the sake of comparison, we adopt the values used in [I], namely 





f h cti 


LIGO 


150 Hz 1.7 x 10~ 22 


Advanced LIGO 


150 Hz 1.3 x 10~ 23 


LISA 


3.88 10~ 3 Hz 1.8 x 10~ 22 



In the figures below we plot the amplitude h(f c ) and the "SNR=1" level . Dividing the first 
by the latter gives an estimation of the SNR. 

3.4.2 Cusp and kink bursts in frequency band of LIGO 

We are now in a position to discuss our predictions for the observed GW burst amplitude 
h(f, N), of both kinks and cusps, as a function of the observation rate N and in the frequency 
band of LIGO, i.e. / = f c = 150Hz. The burst amplitude h(f c , N) is given in ( QUI) together 
with (1251) and ( 126p . It depends on the parameters Gp and p and, in the case of kinks, also on 
the combination qk'. This combination is the fraction q of loops containing junctions multiplied 
by the average number k' of sharp kinks on loops with junctions. As mentioned above, we are 
particularly interested in models for which qk' ^> 1. Therefore to quantify the effect of kink 
proliferation due to the presence of junctions on the GW burst signal we plot, in Fig El Ziasa 
function of N for a range of different qk', with fixed values of the remaining parameters Gp and 
p. We also show h(f c ,N) for cusp bursts, which are not affected by the presence of junctions 
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in our model, as well as the sensitivity limits of LIGO and Advanced LIGO. A discussion of 
how observations of this kind might enable one to estimate the parameter values is given in the 
next subsection. 

As seen in Fig. [2J the proliferation of kinks on loops with junctions significantly increases 
the amplitude of the GW burst signal of kinks and therefore the SNR level of bursts of this 
kind. Even for modest values of qk' ~ 100, the presence of junctions makes the GW bursts 
from kinks observable with LIGO when Gfi> 10 -9 (see also Fig [3] which shows the strongest 
bursts as a function of Gfi). For larger values of qk' one should not only be able to observe the 
largest amplitude bursts from kinks, but also more distant bursts with lower amplitudes (and 
higher rates) hereby obtaining a curve of h as a function of N that can be compared with h(N) 
obtained from theory. 



G\i = 2. 10~ 9 , a = T GyL, T = 50, /> = 0.01 
-20i 




Figure 2: Amplitude of the GW burst signal as a function of the rate N in the frequency 
band of LIGO/Advanced LIGO (f c = 150Hz), for p = 0.01 and G/i = 2 x 10~ 9 . For cusps on 
standard loops (red dashed line) we have taken c = 1; with these parameters z eq ^> z™ sp ;§> 1 
and from (1231) the slope of the cusp curve is —8/11. For kinks (solid black lines) the different 
curves correspond, from bottom to top, to qk 1 = 1, 10 2 , 10 4 , 10 6 , 10 8 and 10 10 . For the large 
values of qk', z^ mk <C 1 so that from fl23|) the slope of the curves is —1/3. For qk' = 1, 10 and at 
large rates one also begins to enter the regime in which z eq ^> ^ nk ^> 1 and the slope changes 
to —7/10. The horizontal solid blue lines correspond to the sensitivity (SNR = 1) of LIGO 
(upper curve) and Advanced LIGO (lower curve). The horizontal dashed (grey) lines denote 
the maximum amplitude /i£nk °f the kink bursts, for q = 1 (upper curve) and q = 10~ 3 (lower 
curve). The vertical line that would indicate the confusion noise limit discussed in Section I3~o1 
lies outside the range in this figure as it corresponds to a rate N = 150Hz=10 9 yr _1 . 

As discussed earlier however the lower limit z c on the redshift, needed to ensure the validity 
of a statistical approach to interpret observations, gives rise to a /c'-independent maximum 
amplitude h ^^ . This is indicated by the horizontal dotted lines in Fig El for two different 
values of q (the upper line has q — 1, the lower one q = 10~ 3 ). For the values of Gfi and p 
taken in Fig [21 one sees that for q — 1, is comparable to the amplitude of the strongest 
(detectable) bursts from cusps and is reached for k' ~ 10 7 . 

The maximal amplitude of bursts reaching LIGO in a reasonable observation time (say with 
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a rate larger than one burst per year) as a function of the tension is plotted in Figure El for 
p = 1 (left panel) and p = 0.01 (right panel). Since the strongest bursts come from the smallest 
accessible redshift and have the smallest accessible rate (see Fig [5]), for cusps with c = 1 (red 
curve), the maximal amplitude is then given by (I25j) with / = f c = 150 Hz and N = 1 yr" 1 . 
For kinks, the maximal amplitude is achieved either for N = 1 yr _1 when < 1 yr" 1 , 

i.e. when Gfi is large enough and given by (1261) (solid black lines); or for N = Ngj^ when 
^kink > 1 J 1 ^ 1 f° r smaller tensions. In the latter case the maximal amplitude is given by (ESI), 
which is represented in Fig [3] by the grey dashed lines in each panel, with q = 1 for the top 
curve and q = 10~ 3 for the bottom one. 

For small values of G/z<10 -10 , FigElshows that for all values of k' cusps remain the strongest 
source of GW bursts that are detectable with LIGO, independently of the reconnection proba- 
bility p and the fraction of loops with junctions q (recall that q < 1). The presence of junctions, 
therefore, does not widen the range of string tensions for which cosmic strings are detectable 
through a GW burst of any type. Also, even though for G[i > 10~ 9 , /i^nk can exceed the ampli- 
tude of the strongest bursts from cusps (for all reasonable p but provided q is not too small), 
we will see in Section H] that in this regime of tensions pulsar timing observations place a rather 
stringent upper bound on kl that essentially rules out those values for which this happens. 




Figure 3: The maximum amplitude of bursts seen by a typical observer at a random location in 
the LIGO/Adv LIGO frequency band during an observation time of one year, as a function of 
Gp,, for p = 1 (left panel) and p = 0.01 (right panel). The red dashed curve gives the amplitude 
of cusps, with c = 1. The amplitude of kink bursts is indicated by the solid black curves for 
(from bottom to top) qk! = 1, 10 2 , 10 4 , 10 6 and 10 8 . The kink curves have an endpoint due to 
the lower cutoff on the redshift. For smaller tensions and the same value of qk', the amplitude 
of kink bursts is given by the grey dashed curve, which only depends on q (with q = 1 for the 
left curve and q = 10~ 3 for the right one). In this case the observation rate N is larger than 
1 yr _1 as discussed in the text. Finally, the horizontal blue lines correspond to the sensitivity 
(SNR = 1) of LIGO (upper curve) and Advanced LIGO (lower curve). 



3.4.3 Cusp and kink bursts in frequency band of LISA 

We now discuss the prospects for observation of the GW burst signal in the frequency band 
f = f c = 3.88 10" 3 Hz of LISA. 
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Figure 4: Amplitude of the GW burst signal as a function of the rate N in the frequency band 
of LISA (f c = 3.88 10~ 3 Hz). The parameters and the meaning of the different curves are the 
same as in Fig [2] except that here qk' = 1, 10 2 , 10 4 , 10 6 , 10 8 , 10 10 . The slope of the cusp curve 
changes from —1/3 to —8/11; while those of the kink curves is —1/3 changing, for small qk' 
and large rates, to —7/10. The horizontal blue line is the "SNR=1" level of LISA given in 
section 13.4. 1[ One sees kink proliferation on loops with junctions can lead to an increase of the 
SNR by several orders of magnitude, before one reaches the cutoff value that arises from the 
lower redshift constraint z c . The vertical dashed dotted line marks the confusion noise limit 
N = 3.3 x 10~ 3 Hz = 10 5 yr" 1 described in section 13.61 beyond which bursts overlap in the 
detector. 



Fig H] is the analogue of Fig [2] above. It shows the observed GW burst amplitude h(f c ,N) 
given in Eq. (1291) as a function of the observation rate N, for both kinks and cusps, and for 
a range of different values of qk'. Apart from the value of the frequency, the parameters are 
identical to those of Fig El The fact that the burst amplitude is inversely proportional to a 
power of the frequency (see Eqs. (I25p-(r2"6"l)) means the amplitude of the signal in the LISA 
band is much larger compared to the LIGO band. This leads to a much larger SNR, so that 
observations of GW bursts in the LISA band open up a new window in parameter space in 
which cosmic string networks are detectable through their gravitational wave burst signal. As 
before, the lower redshift limit z c leads to a maximum amplitude /iJ2nk which, for G[i ~ 10 -9 , 
is reached for k' ~ 10 5 . As can be seen in Fig HI even larger values of k' lead to a wide range 
of high burst rates in which kinks provide the dominant contribution to the GW burst signal 
provided q ~ 1. 

Fig El in turn is the analogue of Fig EJ it shows the maximum amplitude of bursts reaching 
the interferometer during a reasonable observation time, at a rate N > 1 yr" 1 , as a func- 
tion of Gn and for a range of different qk' (note though that here the kink curves are for 
qk' = 1, 10, 10 2 , . . . , 10 6 ). One sees that kinks on standard loops without junctions, namely 
qk' = 1, should be detectable with LISA for tensions > 10" 12 . The presence of junctions does 
not significantly broaden this range. However, in the region where kinks are detectable, the 
proliferation of kinks on loops with junctions leads almost certainly to an increase of the SNR by 
several orders of magnitude, before reaching its cutoff value that arises from the lower redshift 
constraint z c . We also note that the cusp curve in FigOhas an endpoint, at the (p-independent) 
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value G/j ~ [c _1 to(^o/) 2 ^ 3 (ly r_1 )]~ 3 ^ 5 determined from ( I30p . This has a similar origin as the 
endpoints of the kink curves discussed above: As G\i decreases, z m for cusps goes below z c since 
z m ~ (Gfi) 8 / 9 whereas z c ~ (G/i) 1 / 3 , hereby cutting off the observed curve. 




Figure 5: The maximum amplitude of bursts seen by a typical observer at a random location 
in the LISA band during an observation time of one year, as a function of Gfi, for p = 1 (left 
panel) and p = 0.01 (right panel). The red dashed curve gives the amplitude of cusps, with 
c = 1. The amplitude of kink bursts is indicated by the solid black curves for (from bottom to 
top) qk' = 1, 10, 10 2 , 10 3 , 10 4 , 10 5 and 10 6 . Both the cusp and the kink curves have an endpoint 
due to the lower cutoff on the redshift. For smaller tensions and the same value of qk', the 
amplitude of bursts is bounded by the grey dashed curve, which only depends on q (with q = 1 
for the left curve and q = 10 -3 for the right one). In this case the observation rate N is larger 
than 1/year as discussed in the text. Finally, the horizontal blue line corresponds to the LISA 
sensitivity (SNR = 1). 



The large burst amplitudes generally found in the LISA band means there is a priori a rather 
wide regime in which observations might enable one to determine the values of some or even 
all of the four parameters Gfi, p, q and k' in our problem. (We will see in Section 4, however, 
that the stochastic background constrains the window in which this is actually possible.) To 
estimate the parameters one ought to combine the informational from the amplitudes of the 
observed kink and cusp bursts with the additional constraint (s) obtained from one (or two) 
endpoints of the h(N) curves or from changes in the slope of h(N). The latter approach is 
possible only for Gji > 10~ 9 (see FigH]). For G\i < 10~ 10 we find the observable bursts are all 
emitted at small redshift, so that the second derivative of h(N) essentially vanishes, yielding 
no information. However, in this regime the coordinates of the endpoint (s) of the kink and/or 
cusps curves in the (N, /i)-plane are observable, at least for some values of the parameters, 
which provides one or two additional relation(s) between the parameters (on top of the overall 
amplitude of the curves) . 

8 We make no use of information on the position on the sky of the sources of the detected bursts [57] . 
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3.5 Other bursts 



So far we have only discussed bursts generated by kinks and cusps. We now comment on bursts 
produced at kink-kink encounters (k-k bursts) and when kinks cross junctions (k-j). 

The waveforms of these were calculated in [2] and are given in (El) • For a given frequency and 
redshift these bursts have a lower amplitude than those emitted by cusps and kinks. However, 
this could be compensated by the fact that they are emitted in all directions in space. Let us 
start with kink-kink bursts. Following the same derivation as in subsection 13.31 

dN k - k (f,z) = (l + ^VU^W (37) 

where 

~ qk'2 (p^a-VviC*)- 4 ) ■ (38) 

In analogy with 0331) . one can show that the maximum cutoff on the amplitude of k-k bursts 
on loops with junctions is given by 

U™* = Gfi a- 1/3 {ft )~ V V y /3 - (39) 

For LIGO this is below the instrument sensitivity in all parameter space (the requirement for 
individual kink-kink bursts to be observable is (Gfi) 2 q/p > 10~ 5 , 10~ 8 for LIGO, Advanced 
LIGO). For LISA the condition is {Gfx) 2 q/p > 10~ 18 , so that kink-kink bursts should be ob- 
servable provided G[i and q are not too small. The corresponding cutoff on the k-k rate is 

N™1 = fc'Vo)- 1 (40) 

which can be very large since it is proportional to k' 2 . In fact, when /c' 2 (G/i) _1 > 10 21 (resp. 
10 16 ) in the LIGO (resp. LISA) band, this will be close to or above the confusion noise limit 
discussed below in section 13.61 This implies that individual kink-kink bursts should not be 
observable: the sensitivity of LIGO is not sufficient, and whereas the sensitivity of LISA allows 
for their observation, the bursts are in the confusion noise region. However, as we show in 
Section HJ the superposition of many bursts of this kind leads to a significant enhancement of 
the GW stochastic background. 

It is more difficult to estimate the rate of kink-junction events but because kink-junction 
interactions are far less frequent than kink-kink interactions, one expects the observed rate of 
k-j bursts to be negligibly small compared to the rate of k-k bursts. Therefore, since both types 
of bursts have the same waveform, it seems plausible that k-j bursts are irrelevant from an 
observational point of view. 

3.6 From individual bursts to a stochastic background 

Bursts that are produced at large redshifts have small amplitudes but large rates. Therefore 
they can overlap at the detector where their superposition shows up as a Gaussian stochastic 
background. In [I] it was shown that the bursts contributing to the stochastic background seen 
by an instrument observing at a typical frequency f c are those for which the rate satisfies 

AT > f e . (41) 

This bound arises from the fact that, given the burst waveforms, the duration r of a burst in 
the time domain is of the order of their lower frequency cutoff, imposed by 0(1 — 9 m (f, L, z)) 
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[7]. However, since interferometers have larger lower cutoff frequencies^!, the effective duration 
of bursts in the detectors is in fact r m f~ l . During an interval of r, the number of bursts that 
arrive at the detector is of the order of Nr. Therefore one expects to observe a superposition 
of bursts instead of individual bursts when N > f c . 

In our plots of the burst amplitude versus the rate, the inequality (1411) translates into a 
vertical line. Individually detectable bursts must lie to the left of these vertical lines^. In 
the frequency band of LISA the limit is shown in Fig HI where the confusion noise area is 
given by N > 10 5 /year. In the frequency band of LIGO/Advanced LIGO, the limiting rate is 
N > 10 9 /year so that the vertical line lies outside the range of N in Fig [21 The next section is 
devoted to the computation of the characteristic amplitude of the stochastic background. 



4 Superposition of bursts 

4.1 Characteristic amplitude of stochastic background 

Bursts for which the rate satisfies (1411) overlap at the detector and show up as a stochastic 
background rather than as individual bursts. A stochastic background of gravitational waves 
is conventionally characterized by the spectral energy density of gravitational waves, 

"-M - vM (42) 

where p gw is the energy density in GWs and p c the critical density of the universe. For our 
discussion it will be useful to introduce the (dimensionless) characteristic (strain) amplitude of 
the background, h c (f), defined as 

(hVW)) = f S(f - f) h 2 c (f) (43) 

where, as in the remainder of this paper, h(f) is the (dimensionless) logarithmic FT of the 
strain. In terms of h c (f), one finds Q gw (f) is given by 

*W/) = ^(ftofhKf). (44) 

Following prj, the characteristic amplitude h c (f) of the background generated by the super- 
position of bursts is obtained by integrating the square amplitude of individual bursts emitted 
at redshift z, weighted by the number of bursts emitted at that redshift: 

Jmax(z c ,z b ^ b ) J Z 

where N(f, z) is given in ( 1221) and ( 1231) for cusps and kinks respectively. The lower bound on 
the integral is Zb^b, where the subscript stands for 'bursts to background', which is the minimal 
redshift for which bursts arrive superimposed. It follows from ( 1411) that 2r&_>.& is the solution of 

N(zb^b) = f, (46) 



9 An order of magnitude estimate of the cutoff that each instrument applies to different kinds of bursts is 
given by f c , defined in section 13.4. II 

10 As we discuss in Section 2J the existence of a stochastic background has implications for the detectability of 
individual bursts even to the left of this line, since it acts as a "self confusion" noise that could hide the bursts. 
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provided of course that this is larger than z c . Note that for bursts emitted by kinks, z^b 
depends on /, a, p, qk'. 

The integral f!45p must be restricted to values of z satisfying 9 m (z) < 1, where 9 m is given 
in (j7j), to ensure it represents the sum of high frequency bursts. Since 9 m (z) is an increasing 
function of z, this leads to an upper (high frequency) bound Zhf for the integral given by 

9 m (z hf ) = 1. (47) 

By inserting 0161) in (JXJ) and then inverting (JTj), one easily obtains the smooth interpolating 
approximation 

z hf (a, /) = 9 (/ - (1 + ft a) 2 (l + ^ . (48) 

The Heaviside function implies that Zhf — > at f cu t = Hence, for a given value of a all 
spectra are cut at f cut regardless of the type of burst involved. 

Until recently, only the cusp background had been calculated under the assumption that 
this was the leading background. Equation 4.8 in [2Sj gives an analytic approximation to h™ sp 
valid for small frequencies. With our choice of parameters this reads 

h c c nsp « 10(G/i)«- 2/ V 1/2 (^o)" ?/6 - (49) 

Recently, however, it was shown [B] that the amplitude of the kink background is comparable 
to that of the cusp one, even with the standard assumption that the number of kinks per loop 
is of order 1. This calculation relies on a numerical computation of the integral appearing in 
( 145 ft as well as on an analytical approximation of both contributions. 

In the next subsection, we approximate analytically for all frequencies the kink and kink- 
kink backgrounds and we show that even for modest values of k', the latter gives the dominant 
contribution to the GW background. 

4.2 Different contributions to GW background 
4.2.1 The kink contribution 

For kink bursts, using eqs © and ( 1231) . the characteristic strain ( l4"5"j) reduces to 



-1/2^-5/6,, f ,-4 /3( , M/H ( l ZKS ( 1 + ^) 4/3 



1/2 



\Jmax(z^,z^) (1 + Zf^ ) 

We show in Appendix [B] that in the regime in which Zhf ^> z 1 the integral is not sensitive to 
variations of the lower bounds^ and scales as zf/^/z^ 3 . Using (I48p . we see that the condition 
Zhf ^> z 2 eq translates into at f 3> z^ 2 . This yields 

3/2 

h *nk w 1Q ,-V2 ( ^ )a -l/2 ?) - 1 /2 ( ^ ) l/2 (to/) - 1 for / > ^ (51) 



11 One should stress that this is very different from the cusp case studied in [J. One of the important 
observations of that work was that the calculation of h c had to be done by integrating from z^j, (instead of 
from z = 0) in order to avoid counting non overlaping bursts as part of the stochastic background. Though 
this remains obviously true for the kink (and kink-kink) background, our study shows that integrating from 
z = would only introduce a negligible error in the computation of h^ k k . This is because at large frequencies, 
j^k/k k . g (jojjjjj^ted by bursts produced around Zh / and at low frequencies by bursts produced around z = 1 so 
the exact value of Zb_>b does not matter. By contrast, in the case of cusps, starting the integral at z = docs 
introduce an important error in certain regimes. 
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Inserting this into ( PHI) shows that the spectrum is flat at large frequencies and the amplitude 
scales as G^/p (see also equation (63) of [6]). 

We now need to study the regime in which 1 <C at f <C z^ 2 or, in terms of z h f, 1 <C z h j <C 
z 2 q . The key point here is that both z kmk and z k ™ k will be in the z < 1 region. This was 
already discussed in section I3~31 for z kmk (and is true provided that q is not too small). Using 
(USD and ([23]), we see that the condition < 1 translates into (at /) 4/3 < lO^g^a" 1 ?)- 1 
and since we are interested in at f <C z 3 / 2 , a necessary condition for z k ™ k to be 1 is 
10 6 ~ 10~ 2 z 2 q <C (qk')a~ 1 p~ l which is satisfied for all values of the parameters. We are 
therefore in the situation described in the appendix where the main contribution to the integral 
comes from z ~ 1 and the value of the integral is independent of the bounds. This then trivially 
leads to 

= 10(Gp)p^ 2 a- 5 /\t f)-^( q k') 1/2 for — — < f < (52) 

It will be useful to write an interpolating function approximating h kmk at all frequencies 

h^ k « 10(Gp)a 1/2 p~ 1/2 (qk') 1/2 Q (/ - (1 + at,f)~^ (l + ^ ^ . (53) 

An important observation is that h kmk is proportional to k' 1 ! 2 so that, unlike in the case 
of bursts, arbitrarily long kink proliferation would make the amplitude of the background 
arbitrarily large. 



4.2.2 The kink-kink (kk) contribution 

Kink-kink bursts have been discussed in section 13.51 In this case h c is given by 

h k ~ k = vW(G/,) a -v i/2 (*o/)- 3/2 (^ 2 ) i/2 [ r k , {^^] /2 dz 

\Jm&x(z c ,z*J) \ 1 + Z ) 

The derivation of an analytic expression for h k ~ k proceeds in a similar way as the kink case 
above. 

Following the discussion in Appendix [Bj we first note that in the regime in which z^f 3> z^J 2 
the integral is independent of the value of the lower bound and equal to Zhf/z^ 2 = z'^at^f). 
This leads to 

h k - k « 10z- 1 / 2 (G / i)«^Vi/ 2(g ^)i/ 2(to/) -i for / > (55) 

As in the kink case, Q g f is flat at large frequencies. When 1 -C at^f <C z eq , both z k ~ k (which 
is equal to z kmk because it only depends on the population of loops with junctions) and z k ^ 
are < 1 so that 0551) is dominated by z ~ 1 and is of order 1. This yields 

h\- k « 10(G/i)a 1 /V 1 / 2 (gA;' 2 ) 1 / 2 («t /)" 3 / 2 for — < / < ^ . (56) 

ato ato 

Once again, we can summarize our results using a smooth interpolating function 

h k ' k « 10(Gp)a 1/2 p~ 1/2 (qk' 2 ) 1/2 Q (/ - -1) (1 + at /)~ 3 / 2 (l + ^ ^ . (57) 

Note that increasing g/c' 2 simply increases the amplitude of the spectrum whilst leaving its 
frequency dependence unchanged. Decreasing a means decreasing the amplitude while shifting 
the spectrum to higher frequencies. 
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4.2.3 Comparison with cusp background 

Dividing (JF3) by §2]), one obtains h k ~ k /h kink = y/k>$(at Q f) where = (1 + x) 1/6 (l + 

x / z eq )~ l l 2 (l + x/zf/ 2 ) 1 / 3 . This function is of order 1 for all values of x so we conclude that 
h k ~ k is always larger than h k c by a factor of order \fk' . Note that this also implies that in the 
absence of kink proliferation, when q = k! = 1, kink-kink events contribute at the same order as 
kink bursts to the stochastic background. 

Furthermore, it was shown in [5] that the contribution from cusps to the stochastic back- 
ground is of comparable order to the standard kink one (q — k' = 1). We therefore conclude 
that as soon as qk! ^> 1, the dominant contribution to the GW background comes from the k-k 
bursts. For this reason, from now on we mainly focus on the k-k background. 



4.3 Self-confusion noise 

We now compare the amplitude of the bursts with that of the dominant k-k contribution to the 
background in order to understand to what extent the latter acts as (self) confusion noise for 
the detection of individual bursts. In Fig E] we plot both the amplitude of the k-k background 
for a range of different values of qk' 2 together with the maximal amplitudes of bursts, as a 
function of Gfi and in the frequency bands of LIGO and LISA. The maximal amplitude for the 
bursts is taken to be the g-dependent cutoff amplitude given in Eq. fj33j) for the kink bursts, 
and the amplitude corresponding to a rate of 1 event per year for cusps. 



a = r<7|i, r = 50, p=\, 

f=fLIGO=™ Hz 



a = T G\i, T = 50, p=\, 

f C =f L ISA = ^ 10 
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Figure 6: The solid red lines give the amplitude of the stochastic gravitational wave background 
generated from the superposition of kink-kink bursts as a function of G/i, in the frequency band 
of LIGO (left) and LISA (right) for different values of qk' 2 (from bottom to top: qk' 2 = 1, 10 3 , 10 6 
and also 10 9 in the LIGO band). For comparison, we also plot the maximum amplitude of 
the individual bursts. The two parallel black dashed lines in both panels denote the cutoff 
amplitudes for kink bursts (for q — 1 and q = 10~ 3 ) and the long-dashed black line is the 
amplitude of cusp bursts arriving at the instrument at the rate of 1/year. One sees that for 
large but reasonable values of qk' 2 , the GW background hides the individual high frequency 
bursts. 

Remarkably, the values of qk' 2 required for the background to dominate over the bursts are 
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not particularly large, typically around 10 6 . Therefore there is a parameter regime in which 
kink proliferation on loops with junctions makes individual bursts simply unobservable. 



5 Observational constraints on GW background 

We end with a discussion of the constraints on the different parameters of our cosmic-string 
model arising from observations of millisecond pulsars and LIGO, as well as the future prospects 
of observation with LISA. After a brief review of the current observations we summarize the 
resulting constraints in Section 15.41 and discuss their implications for the GW burst signal in 
our model. 



5.1 Milisecond Pulsars 

The latest constraints on the dimensionless strain of gravitational waves imposed by millisecond 
pulsars are derived in [29J. For a frequency dependence of the form 

h 6 (f)=A(J-^ ' (58) 

in the domain [1/20 yr _1 , 1 yr _1 ], the dimensionless amplitude A is constrained to be 

A<10~ 15 . (59) 

On comparison with the string spectrum calculated in (15"7|) . it follows that ( 1351) can be used to 
apply a constraint provided the frequency interval [(ato) -1 , z eq (ato)~ 1 ] intersects yr _1 , 1 yr -1 ]. 
In practice, for the range of values of a we consider, this reduces tern (ato)^ 1 < 1 yr _1 
i.e. roughly a > 10 -10 . 

Under this assumption and using Eqs. (j5"7j) - (j5"91) . we obtain the constraint p~ 1 ^ 2 (qk' 2 ) 1 ^ 2 < 
10~ 16 (1 yr)- 3 / 2 t^ /2 r which, since p < 1 and taking T = 50, yields the final order of magnitude 
constraint coming from pulsars 

qk' 2 < 25 for a > 10" 10 i.e. Gfi > 10" 12 (60) 

Note that within the range a > 10 -10 , the constraint on parameter space is independent of a. 



5.2 LIGO 

The ground based interferometer LIGO has already begun to take data. An upper limit on the 
amplitude of a stochastic background of gravitational waves based on the data from a two year 
science run, assuming a frequency dependence in the frequency band of LIGO (~ 10 2 Hz) of the 
form ^ 

*W/) = fi 7 (io/^J around 102 Hz ( 61 ) 

is given in [HD]- For the k-k background and for the range of values of a of interest, 7 = 0. 
Indeed, h c oc f- 1 for / > fa an d if a > 10" 16 , fa < 10 2 Hz so Q gw oc f 2 h 2 c (f) oc p with 7 = 
in the frequency band of LIGO. 

12 Thc other condition, namely z eq (ato)^ 1 > 1/20 yr _1 only imposes that a < 10 -5 . We also note that the 
slope of the spectrum can never be —1 instead of —3/2 in this frequency interval for this would require values 
of a larger than 10~ 5 . 



20 



The 95% confidence upper limit on fl gw that applies to the k-k background reads 

n gw (f) <tt = 6.9 1(T 6 around 10 2 Hz (62) 
which, using (jHJ) and H = 72 km/s/Mpc translates into 

h c (f) < 2.4 KT 21 f-y-J around 10 2 Hz . (63) 

This yields, in our parameter space, the constraint p~ x qk' 2 a < (2.4 1CT 21 ) 2 (^) t^z eq , i.e. using 
T = 50 and p < 1, 

qk' 2 a < 0.1 (64) 

valid for 

a > 10~ 16 . (65) 

Note that the constraint on fl gw ( l6"2"j) is expected to be improved by 4 orders of magnitude 
with Advanced LIGO, which means an improvement by 2 orders of magnitude on the constraint 
on h c fl6"3"|) and an improvement by 4 orders of magnitude for the constraint on aqk' 2 . 



5.3 LISA 



In the not-too-distant future, hopefully the space interferometer LISA will enable us to probe a 
different frequency band of gravitational waves around ~ 10~ 3 Hz. Using the sensitivity curves 
for the strain h c integrated over a year provided in [21] with the default values of the parameters 
of LISA, one finds LISA should be able to detect the k-k background for a as small as 10~ 14 
even in models where qk' 2 is not much larger than 1, as illustrated in Fig 15.31 
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Figure 7: The future space interferometer LISA should be able to detect the GW background 
generated from bursts emitted at kink-kink encounters for string tensions as low as 10 -15 and 
for all values of qk' 2 > 1. 
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5.4 Discussion 



The observations of millisecond pulsars and LIGO place rather stringent constraints on the 
different parameters of our cosmic-string model which are summarized in Fig 15.4} where we 
plot the various constraints in the (Gfi, qk' 2 ) plane. 
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Figure 8: Summary of the different observational constraints on the parameters Gfi and qk' 2 . 
For tensions Gfi > 10~ 12 pulsar observations provide the most stringent constraints and enforce 
qk' 2 to be less than ~ 10 2 . For smaller values of Gfi, LIGO places an upper bound on qk' 2 but 
this is very high. With Advanced LIGO one will be able to significantly improve on this. LISA 
in turn should be sensitive to all value of (Gfi, qk' 2 ) in the range given here. 

These observational constraints are relevant for our discussion in Section 13.41 of the GW 
burst signal in the class of models under consideration. Given the above assumptions about 
the network evolution and also ([HD and ffTTj) . we believe the most important astrophysical 
implications of this analysis are the following: 

• For string tensions Gfi > 10 -12 , the pulsar data imply qk' 2 < 10 2 . Hence models with an 
important kink proliferation on a significant proportion of the loops in the network are 
ruled out in this range of tensions. (Note that a similar conclusion applies to a standard 
network of cosmic strings with no junctions and hence no kink proliferation. If the tension 
of these strings Gfi > 10~ 12 then the number of kinks k on the loops must satisfy k < 10.) 

As a consequence, in the remaining models in this range of tensions, the predictions of 
[I] for the gravitational wave signal essentially apply. In particular, cusp bursts remain 
the dominant source of high frequency gravitational waves. Kink bursts should not be 
observable with LIGO or Advanced LIGO except for large tensions and provided p is 
small. With LISA, bursts from kinks should be observable but with an SNR that is much 
smaller than the SNR of cusp bursts. The observation of kink bursts however should help 
to estimate the parameters values. 

• For string tensions Gfi < 10 -12 , the current observational constraints on qk' 2 are much 
weaker and models exhibiting significant kink proliferation (qk' 2 ^> 1) are allowed. If this 
is the case, the consequences on the GW predictions are the following. 

The SNR of kink bursts increases but only very slightly (never by more than one order of 
magnitude, see Figs. ElandE]), because one quickly reaches the cutoff value given in (13111) . 
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In particular, individual kink bursts remain below the sensitivity of LIGO and Advanced 
LIGO. They will be accessible to LISA, with low SNR, only if G/i is close to 10~ 12 and p 
is small. 

Furthermore, for large values of qk' 2 the GW signal from individual bursts in the LISA 
band is hidden by the confusion noise generated by smaller bursts from kink-kink encoun- 
ters. Indeed, for values qk' 2 > 10 4 , the GW background generated from the superposition 
of kink-kink bursts becomes so strong that it hides the individual sharp bursts from cusps 
and kinks (see Fig El) . 

Hence in this parameter regime the most important astrophysical implication of the pres- 
ence of junctions is the absence of observable individual bursts and the presence of a 
strong GW background detectable by LISA and maybe by LIGO, with a spectrum given 
in (JSZI). 

6 Conclusion 

A network of cosmic string loops generates a strongly non- Gaussian GW signal which includes 
sharp bursts at cusps and kinks that stand above the confusion noise generated by many smaller 
overlapping bursts. For loops without junctions, the bursts at cusps provide the dominant con- 
tribution to the GW signal in all frequency bands of observational interest. However, in the 
presence of junctions the number of sharp kinks grows exponentially, leading to an average 
number k! of large amplitude kinks on loops with junctions that can be several orders of mag- 
nitude larger than that is expected on loops without junctions. In this paper we have studied 
the effect of this proliferation of kinks on cosmic string loops with junctions on the GW burst 
signal emitted from string networks of this kind. 

We have calculated the rate of occurrence and the distribution in amplitude of the GW 
bursts emitted at cusps and kinks in the frequency bands of LIGO and LISA, as a function of 
the string tension Gfi, the number k' of sharp kinks on loops with junctions, the fraction q of 
loops that have junctions and the reconnection probability p. To do so we have used a simple 
'scaling' model for the cosmic string network. Our predictions for GW observations depend on 
the simplifying assumptions involved in this model. These include, in particular, that loops are 
characterized by a single length scale L (see Eq. flS])), and that all loops formed at time t have 
the same size L ~ at, with a = 50G/i. 

The amplitude of individual kink bursts depends on the combination qk'. In models with 
junctions one generally expects qk' ^> 1, whereas in models without junctions one usually 
assumes loops have only one or a few kinks. We have seen that this difference lowers the redshifts 
required to achieve observable rates and therefore can lead to a significant enhancement of the 
GW burst signal of kinks in loops with junctions relative to standard loops for given values of 
the string tension and p. 

In fact, one might have thought that by taking qk' sufficiently large, one could easily make 
the SNR level of kink bursts much larger than the SNR of cusps. However, we have seen there 
is an important limitation that arises from the density of loops in the model. This is because 
below a certain redshift the probability of finding a loop - and therefore a GW burst- becomes 
exceedingly small. Since the largest amplitude bursts are emitted at the smallest redshift this 
translates, in the scaling model, in an upper bound on the amplitude of kink bursts. The latter 
depends on G/i, p and q but lies, for string tensions below 10~ 10 , always below the amplitude 
of cusp bursts both in the LIGO and LISA frequency bands. For larger string tensions this 
is a priori not the case. However, in this case significant kink proliferation is ruled out by 



23 



millisecond pulsar observations, so that cusps remain the strongest source of gravitational wave 
bursts in both bands. 

Another effect of kink proliferation is to significantly enhance the stochastic background 
generated by the superposition of overlapping bursts at the detector. In Section [H we computed 
analytically the spectrum of the kink and of the kink-kink backgrounds. We showed that as 
soon as qk' 2 > 1, the latter becomes dominant over the cusp and the kink backgrounds. 

By combining our results on the stochastic background with the latest observations from 
millisecond pulsars and LIGO, we were able to constrain the parameters in the model. For 
string tensions G[i > 10~ 12 , we find significant kink proliferation is not allowed and hence as 
far as the GW predictions are concerned, the results of [I] essentially apply. For G\i < 1CT 12 the 
observational constraints on kink proliferation are much weaker. In models with a large number 
of kinks on loops, the SNR of individual kink bursts in LISA or LIGO increases though only 
by one order of magnitude at most because of the cutoff mentioned above. Most importantly, 
however, the stochastic background generated from the incoherent superposition of kink-kink 
bursts is strongly enhanced and becomes dominant over the cusp and kink backgrounds when 
qk' 2 > 1, and it even hides all individual bursts from both cusps and kinks for qk' 2 > 10 4 . 
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A GW signal dominated by large amplitude kinks 

The goal of this appendix is to illustrate the fact that GWs produced by kinks on loops with 
junctions are dominated by large amplitude kinks (yl>0.25), despite the larger number of 
small amplitude kinks (A <C 1). (Recall that whenever a kink propagates through a junction 
it generates 3 daughter kinks whose amplitude, is generally smaller than that of the initial 
kink [3J.) As a result, when computing the bursts rates or the characteristic amplitude of the 
stochastic background produced by the incoherent superposition of many bursts, we will not 
need to use the full distribution of kink amplitudes on such loops but instead only the number 
k' of kinks that have an amplitude >0.25 as done throughout the paper. 

Let f{A) be the distribution of kink amplitudes on a loop at the end of the proliferation 
phase (so that f(A)dA is the number of kinks that have an amplitude between A and A + dA). 
In principle, f{A) depends on the initial conditions of the particular loop we are considering. 
But, as was argued in [3], the details of these initial conditions are lost after a few time intervals 
of L (the typical size of the loop). We will therefore consider f(A) to be the same for all the 
loops in the network. 

The form of f(A) was computed numerically in (3] for an idealized loop, using a tree con- 
struction. Apart from the string tensions, the main parameter in that setup is the number n 
of generations for which proliferation takes place (physically its time duration in units of L). 
Simulations of realistic loops lead us to consider n of the order of 10 but this parameter could 
be much larger. From now on we will use the notation f n (A) when we want to indicate the 
particular value of n that we used, and simply f(A) elsewhere. 
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The total number of kinks obviously increases (exponentially) with n but f n {A) becomes 
peaked around some small value that goes to zero exponentially fast with n: the vast majority 
of kinks have a very small amplitude. However, the number of kinks that have an amplitude 
larger than some fixed value also increases exponentially. 

A.l Burst rates 

Intuitively, small amplitude kinks produce bursts that can only be observed experimentally if 
they were produced at small enough redshift. The question of their negligeability is therefore a 
competition between their large number and the fact that the number of loops decreases with 
redshift. 

To be more precise, let us determine the rate of events N(h) that have an amplitude larger 
than some fixed value h (that could be for instance the experimental threshold). 

To answer this question, we need to slightly modify formula ([6]) to take into account the fact 
that we no longer consider kinks of amplitude A = 1 exclusively. The amplitude of a GW burst 
produced at a kink of amplitude A on a loop at redshift z (and observed at some frequency /) 
is simply 

h ^z)=A ^ L \ +Z (66) 

((i + z)Lfyi i t z 

where we have now removed the G function that appeared in because for reasonable values 
of / (in the band of LISA or LIGO for instance) it only imposes a higher cutoff value on z that 
is much larger than all redshifts of interest for bursts. We will also need an updated version 
of (TTU|) to compute the rate of bursts. We now define dN(A, z) the number of bursts reaching 
us and that were emitted between z and z + dz from kinks with an amplitude between A and 
A + dA. This quantity reads 

dN(A, z) = + z)- 1 qn(z)f{A)dV(z)dA (67) 

We can simplify the expressions above by remarking that we are only interested in bursts 
coming from the region z -C 1 (see section 3 of this paper). In this region, the density of loops 
n(z), the beaming angle 9 m and the typical size of the loops L do not depend on z and can 
be approximated by their values today. We could express those in terms of our parameters 
p, a but we are only interested here in the form of the dependence on A and z so we will 
simply use the notations no, Oq and Lq. Furthermore, we will use dV(z) ~ t\z 2 dz. Defining 
C = GuL 1/3 /~ 2 / 3 to 1 and D = O^L^qu^ we obtain 

A 

h(A,z)=C- and dN(A,z) = Df(A)z 2 dzdA (68) 

z 

It is then clear that N(h) is obtained by integrating dN(A, z) on the domain T>(h) shown 
below. 
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N(h) 



dN{A,z) 



D / f(A) 
Ja=o 



h 



z 2 dz dA 



DC 3 



f{A)A 3 dA (69) 



IV{h) JA=0 ' \Jz=0 ) 3/i 3 JA=0 

The question is now to know if this last integral is dominated by large values of A or by the 
A <^ 1 region. We define for convenience the funtions 

F k , n (x)= f 1 f n (A)A k dA (70) 

J X 

We see from figure IA~T1 that the dominant contribution to f\ =0 f(A)A 3 dA comes from A « 1 
(as x starts decreasing from 1, the red curve quickly reaches a plateau). In that figure, we used 
n = 13. 




Figure 10: As x decreases from 1 to 0, ^3,13(0;) (red) strongly grows and quickly reaches a 
plateau which shows that the main contribution to the integral appearing in fl6"9"|) (namely 
-^3,13(0)) comes from Ami. By contrast, the main contribution to the kink distribution 
(illustrated by -^0,13) comes from smaller amplitudes. A good approximation to the integral in 
is F 3 , 13 (0) w'Fo,i 3 (1/4) 



Since the dominant contribution comes from A « 1, we can estimate the integral by replacing 
A 3 by 1 and simply counting the number k' of kinks that have an amplitude of the order of 1, 
i.e. between 1 and say 1/4 

f 1 f(A)A 3 dA w f 1 f(A)A 3 dA w l 3 f{A)dA = k' (71) 

JA=0 JA=l/i JA=1/A 
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We see from figure lA~T1 that such an estimation is a very good approximation at least for our 
order of magnitude calculations: the intersection between both curves occurs close to x = 1/4 
and the red one has already reached its plateau value. We checked that this was not very 
sensible to a different choice of n. Of course, the value of k' does depend on the number of 
generations n and we leave it as a free parameter. 

Therefore, only large amplitude bursts contribute to the rate of bursts we might observe 
and the picture that we have used in this paper, namely loops with k' kinks that all have the 
same amplitude A = 1 yields a very good approximation for the rate of bursts. 



A. 2 Stochastic background 

We now address the question of the importance of small amplitude kinks in the background 
produced by the incoherent superposition of many bursts? The (squared) characteristic am- 
plitude of such a background is computed by summing the squared amplitude of bursts that 
actually overlap at the detector. 

In situations (as the one we consider throughout this paper except for this appendix) where 
there is a direct correspondance between the amplitude of the bursts and the redshift (and also 
between rates and redshifts), this sum can be expressed as in ( 1451) as an integral over redshifts 
larger than the value z^b for which the rate becomes larger than the inverse timescale of the 
experiment. If we allow for a whole distribution of kink amplitudes instead of considering 
that all kinks have amplitude 1, the situation becomes more complicated. We use the same 
procedure as [5] to remove the infrequent (non overlapping) large amplitude bursts, namely we 
only include in the sum bursts that have an amplitude h smaller than the value for which 
N{h) > f. This amounts to performing the integral of the quantity f~ 1 h 2 (A, z)dN(A, z) over 
the domain T>' showed in figure IA.2I 



h(A, z) = h* 




^ Z hf ' Z 

Figure 11: Domain of integration of f~ 1 h 2 (A, z)dN(A, z) required to compute h c . 



Let z(A) be defined by h(A, z(A)) = h*. If z(A = 1) C 1 (which is the most frequent case) 
then this is a stright line in the (A,z) plane. The integral can be written 

hi = /' f{A)A 2 ( f h ' d A^pA h * {A = !, z)dz ) dA (72) 

JA=Q \Jz(A) J J 

Except for its bounds (and obviously for a factor k'), the inside integral over redshifts turns out 
to be exactly the one that we encountered for the calculation of h c in HSJ We already showed 
that it was never sensible to the variations of its lower bound. We are therefore simply left with 
Jjj =0 f(A)A 2 dA as an overall multiplicative factor. As in the case described above of bursts, 
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this integral is dominated by A « 1 and an be approximated by k! = f x u f(A) the number of 
kinks of large amplitude. The validity of this approximation is illustrated in figure IA.2I 
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Figure 12: As x decreases from 1 to 0, ^2,13(2?) (red) strongly grows and quickly reaches a 
plateau which shows that the main contribution to j\ =0 f(A)A 2 dA comes from A w 1. A good 
approximation to j\ =0 f(A)A 2 dA is F 0t i 3 (l/A) i.e. the number of large amplitude kinks that we 
studied in [3] 

Here again, simply considering that loops have k! large amplitude kinks instead of consid- 
ering a whole distribution of kink amplitudes yields the correct result for the background. 



B Redshift integrals 

We provide an analytical description of the integrals that appear in the calculation of the 
different types of backgrounds. We will always assume here that z a <C z h . 

B.l The kink-kink integral 

In this appendix, we study the behaviour of the function 



In order to give an idea of where the main contributions to this integral come from, we plot the 
function zfkk(z) in figure lETTl 

This function has a local maximum around z — 1, a local minimum around z = z eq and 
then diverges (as -§j?) when z goes to infinity. The value at which it becomes larger than the 



local maximum is z ~ z'H 2 

Using the fact that for z< 1, fkk{ z ) ~ 1, for 1 <C z <S z eq , fkk{ z ) ~ z~ 3 ^ 2 and for z 3> z eq , 

3/2 

fkk(z) ~ z~' q l 2 it is easy to see that J Q eg f(z)dz = 0(1) and that the dominant contribution 




3/2 



(73) 
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to this integral comes from an interval around z — 1. Therefore, in our two regimes of interest 



z a <C 1 and 1 C Z(, « z^ 2 : then g doesn't depend (at leading order) on z a or z b . 
g(z a ,z b ) « 1 

z a ^ z ei 2 an d ^6 ^ ^eg 2: then J Z y 2 f(z)dz pa ^> 1 dominates over the rest of the 

^eq %eq 

integral and g(z a , z b ) « -f^. 
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Figure 13: In the case of the kink or the kink-kink backgrounds, the integral is either dominated 
by large redshifts or by z ~ 1. The situation is very different for cusps where the domination 
comes from the smallest redshifts (or from z ~ 1) 

B.2 The kink integral 

The integral appearing in the calculation of the kink background is 



i \ f" U + 

9kink — [Za,Zb) ~^ / 7T 

J z„ ( 1 



*»- (1 + z/z e9 ) 



4/3 



(l + ^p73— ^= / fkink(z)dz 



(74) 



As can be seen from figure IB. 1L the shape of zfkink(z) is the same as in the kk case so the 
discussion above applies after correctly changing the exponents. In particular, for our regimes 
of interest 

• z a < 1 and 1 < z b < : g(z a , z b ) « 1 



z a < z^ and z fe > z^: #(z a , z fe ) 



2/3 



eg • 



29 



B.3 The cusp integral 

The integral appearing in the calculation of the cusp background is 

[*> (l+z/z eq ) 7 / 6 m 

9cusp — \z a ,z b ) — > J — z yu/ & J Jcusp\Z)aZ {(0) 

As before, in an effort to identify the main contributions to such an integral, we plot zf cusp (z) 
The situation is very different than in the kk case. There is again a local maximum around 
z — 1 where zf cusp (z) = 1 and then we have zf cusp > 1 for z > z 7 J q 2 . Among the regimes of 
interest, there will be one where both z a and z^ lie in the interval [1,2: J/ 2 ] which explains why 
the value of z a is crucial in the case of cusps. 



References 



T. Damour and A. Vilenkin, "Gravitational wave bursts from cosmic strings," 



Phys. Rev. Lett. 85 (2000) 3761-3764] |arXiv:gr-qc/0004075[ 



[2] 

[3] 
[4] 
[5] 
[6] 

[7] 



P. Binetruy, A. Bohe, T. Hertog, and D. A. Steer, "Gravitational Wave Bursts from 



Cosmic Superstrings with Y- junctions," \Phys. Rev. D80 (2009) 123510 



arXiv: 0907 .4522 [hep-th] 



P. Binetruy, A. Bohe, T. Hertog, and D. A. Steer, "Proliferation of sharp kinks on cosmic 



(super-) string loops with junctions," |arXiv : 1005 . 2426 [hep-th] 



T. Damour and A. Vilenkin, "Gravitational wave bursts from cusps and kinks on cosmic 



strings," \Phys. Rev. D64 (2001) 064008| |arXiv:gr-qc/0104026[ 



L. Lorenz, C. Ringeval, and M. Sakellariadou, "Cosmic string loop distribution on all 
length scales and at any redshift," arXiv: 10 06.0931 [astro-ph . C0"T] 

S. Olmez, V. Mandic, and X. Siemens, "Gravitational- Wave Stochastic Background from 
Kinks and Cusps on Cosmic Strings," \Phys. Rev. D81 (2010) 104028[ 
|arXiv: 1004.0890 [astro-ph. CO] 



X. Siemens et ai, "Gravitational wave bursts from cosmic (super) strings: Quantitative 
analysis and constraints," \Phys. Rev. D73 (2006) 105001| |arXiv:gr-qc/0603115[ 



X. Siemens, V. Mandic, and J. Creighton, "Gravitational wave stochastic background 



[9] 
[10] 

[11] 



from cosmic (super) strings," Phys. Rev. Lett. 98 (2007) 111101} 
|arXiv : astro-ph/0610920| 

J. Polchinski and J. V. Rocha, "Analytic Study of Small Scale Structure on Cosmic 
Strings," \Phys. Rev. D74 (2006) 083504| |arXiv : hep-ph/0606205l 

J. V. Rocha, "Scaling solution for small cosmic string loops," 



Phys. Rev. Lett. 100 (2008) 071601, arXiv: 0709 . 3284 [gr-qc] 



F. Dubath, J. Polchinski, and J. V. Rocha, "Cosmic String Loops, Large and Small," 



Phys. Rev. D77 (2008) 123528, ar Xiv: 071 1.0 994 [a stro-ph] 



30 



[12] 

[13] 
[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[20] 

[21] 
[22] 

[23] 

[24] 
[25] 
[26] 

[27] 
[28] 



A. Avgoustidis and E. P. S. Shellard, "Velocity-Dependent Models for 



Non-Abelian/Entangled String Networks," Phys. Rev. D78 (2008) 103510 



arXiv: 0705 .3395 [astro-ph] 



M. Hindmarsh and P. M. Saffin, "Scaling in a SU(2)/Z(3) model of cosmic superstring 
networks," JEEP 08 (2006) 066, |arXiv : hep-th/0605014| 



A. Rajantie, M. Sakellariadou, and H. Stoica, "Numerical experiments with p F- and q 
D-strings: the formation of (p,q) bound states," \JCAP 0711 (2007) 021, 



|arXiv: 0706 .3662 [hepth] 



S. H. H. Tye, I. Wasserman, and M. Wyman, "Scaling of multi-tension cosmic superstring 



networks," \Phys. Rev. D71 (2005) 103508| |arXiv : astfo~- ph/0503506 



J. Urrestilla and A. Vilenkin, "Evolution of cosmic superstring networks: a numerical 



simulation," JEEP 02 (2008) 037, arX iv: 0712. 1146 [hep-th] 



E. J. Copeland and P. M. Saffin, "On the evolution of cosmic-superstring networks, 
JEEP 11 (2005) 023, |arXiv : hep-th/0505TToj 

T. W. B. Kibble, "EVOLUTION OF A SYSTEM OF COSMIC STRINGS " 



Nucl. Phys. B252 (1985) 227. 



C. J. A. P. Martins and E. P. S. Shellard, "Quantitative String Evolution," 
Phys. Rev. D54 (1996) 2535-2556| |arXiv:hep-ph/9602271| 



C. J. A. P. Martins and E. P. S. Shellard, "Extending the velocity-dependent one-scale 



string evolution model," \Phys. Rev. D65 (2002) 043514| |arXiv : hep-"ph/0003298 



J. Urrestilla, "Private communication,". 

E. J. Copeland, T. W. B. Kibble, and D. A. Steer, "Collisions of strings with Y 



junctions," \Phys. Rev. Lett. 97 (2006) 021602[|arXiv:hep-th7b601153 



H. Firouzjahi, J. Karouby, S. Khosravi, and R. Brandenberger, "Zipping and Unzipping 



of Cosmic String Loops in Collision," \Phys. Rev. D80 (2009) 083508| 
|arXiv : 0907 . 4986 [hep-th] | 

D. Garfin kle and T. Vachaspati, "RAD IATION FROM KINKY, CUSPLESS COSMIC 
LOOPS," \Phys. Rev. D36 (1987) 2229] 

N. Bevis et ai, "Evolution and stability of cosmic string loops with Y- junctions," Phys. 
Rev. D80 (2009) 1 25030, [arXiv: 0904. 2 127 [hep-th]] 

LIGO Scientific Collaboration, B. P. Abbott et al, "First LIGO search for 
gravitational wave bursts from cosmic (super) strings, " \Phys. Rev. D80 (2009) 062002 
|arXiv: 0904. 4718 [astro-ph. CO] [ 

M. I. Cohen, C. Cutler, and M. Vallisneri, "Searches for Cosmic-String 
Gravitational- Wave Bursts in Mock LISA Data," [arXiv: 1002.4153 [gr-qcT 



T. Damour and A. Vilenkin, "Gravitational radiation from cosmic (super) strings: Bursts, 



stochastic background, and observational windows," Phys. Rev. D71 (2005) 063510[ 
|arXiv : hep-th/04102~22~j 



31 



[29] F. A. Jenet et al, "Upper bounds on the low- frequency stochastic gravitational wave 
background from pulsar timing observations: Current limits and future prospects," 



Astrophys. J. 653 (2006) 1571-15761 |arXiv: astro-ph/0609013 



[30] LIGO Scientific Collaboration, B. P. Abbott et al, "An Upper Limit on the Stochastic 
Gravitational- Wave Background of Cosmological Origin, 



Nature 460 (2009) 990 



arXiv: 0910 .5772 [astro-ph . CO]] 



[31] http://www.srl.caltech.edu/ shane/sensitivity. 



32 



